Spin-orbit coupled Fermi liquid theory of ultra-cold magnetic dipolar fermions 
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We investigate Fermi liquid states of the ultra-cold magnetic dipolar Fermi gases in the simplest 
two-component case including both thermodynamic instabilities and collective excitations. The 
magnetic dipolar interaction is invariant under the simultaneous spin-orbit rotation, but not under 
either the spin or the orbit one. Therefore, the corresponding Fermi liquid theory is intrinsically 
spin-orbit coupled. This is a fundamental feature of magnetic dipolar Fermi gases different from 
electric dipolar ones. The Landau interaction matrix is calculated and is diagonalized in terms 
of the spin-orbit coupled partial-wave channels of the total angular momentum J. The leading 
thermodynamic instabilities lie in the channels of ferromagnetism hybridized with the ferronematic 
order with J — 1 + and the spin-current mode with J — 1~, where + and — represent even and 
odd parities, respectively. An exotic propagating collective mode is identified as spin-orbit coupled 
Fermi surface oscillations in which spin distribution on the Fermi surface exhibits a topologically 
nontrivial hedgehog configuration. 



PACS numbers: 03.75.Ss,05.30.Fk,75.80.+q,71.10.Ay 



I. INTRODUCTION 

Recent experimental progress of ultracold electric dipo- 
lar heteronuclear molecules has become a major focus of 
ultracold atom physics 1 - - — . Electric dipole moments are 
essentially classic polarization vectors induced by the ex- 
ternal electric field. When they are aligned along the z 
axis, the electric dipolar interaction becomes anisotropic 
exhibiting the d r .2_ 3z 2-type anisotropy. In Fermi sys- 
tems, this anisotropy has important effects on many- 
body physics including both single-particle and collective 
properties^— . Fermi surfaces of polarized electric dipo- 
lar fermions exhibit quadrupolar distortion elongated 
along the z axis 4 ^" 7 -* 1 ^. Various Fermi surface instabil- 
ities have been investigated including the Pomeranchuk 
type nematic distortions^ 7 - and stripelike ordering a 10 ' 14 . 
The collective excitations of the zero sound mode exhibit 
anisotropic dispersions: The sound velocity is largest if 
the propagation wavevector q is along the z axis, and 
the sound is damped if q lies in the xy planed. Un- 
der the dipolar anisotropy, the phenomenological Lan- 
dau interaction parameters become tridiagonal matrices, 
which are calculated at the Hartree-Fock level^i 7 -, and the 
anisotropic Fermi liquid theory for such systems has been 
systematically studied 7 . 

The magnetic dipolar gases are another type of dipolar 
system. Compared to the extensive research on electric 
dipolar Fermi systems, the study on magnetic dipolar 
ones is a new direction of research. On the experimen- 
tal side, laser cooling and trapping Fermi atoms with 
large magnetic dipole moments (e.g., 161 Dy and 163 Dy 
with /i = 10/is)i£~— have been achieved, which provides 
a new opportunity to study exotic many-body physics 
with magnetic dipolar interactions. There has also been 
a great amount of progress for realizing Bose-Einstein 
condensations of magnetic dipolar atomsi 7 - - — . 

Although the energy scale of the magnetic dipolar in- 
teraction is much weaker than that of the electric one, it 
is conceptually more interesting if magnetic dipoles are 



not aligned by external fields. Magnetic dipole moments 
are proportional to the hyperfine spin up to a Lande fac- 
tor, thus, they are quantum-mechanical operators rather 
than the nonquantized classic vectors as electric dipole 
moments are. Furthermore, there is no need to use ex- 
ternal fields to induce magnetic dipole moments. In fact, 
the unpolarized magnetic dipolar systems are isotropic. 
The dipolar interaction does not conserve spin nor orbit 
angular momentum, but is invariant under simultaneous 
spin-orbit (SO) rotation. This is essentially a spin-orbit 
coupled interaction. Different from the usual spin-orbit 
coupling of electrons in solids, this coupling appears at 
the interaction level but not at the kinetic-energy level. 

The study of many-body physics of magnetic dipolar 
Fermi gases is just at the beginning. For the Fermi liquid 
properties, although magnetic dipolar Fermi gases were 
studied early in Refs. [22[ and [6|, the magnetic dipoles 
are frozen, thus, their behavior is not much different from 
the electric ones. It is the spin-orbit coupled nature that 
distinguishes non-polarized magnetic dipolar Fermi gases 
from polarized electric ones. The study along this line 
was was pioneered by Fregoso and Fradki n 23 i 24 . They 
studied the coupling between ferromagnetic and ferrone- 
matic orders, thus, spin polarization distorts the spher- 
ical Fermi surfaces and leads to a spin-orbit coupling in 
the single-particle spectrum. 

Since Cooper pairing superfluidity is another impor- 
tant aspect of the many-body phase, we also briefly sum- 
marize the current progress in electric and magnetic dipo- 
lar systems. For the single-component electric dipolar 
gases, the simplest possible pairing lies in the p-wave 
channel because s-wave pairing is not allowed by the 
Pauli exclusion principle. The dipolar anisotropy se- 
lects the ^-channel pairing 2 ^— . Interestingly, for the 
two-component case, the dipolar interaction still favors 
the triplet pairing in the p z channel even though the s 
wave is also allowed. It provides a robust mechanism for 
the triplet pairing to the first order in the interaction 
strength 3 ^—. The mixing between the singlet and the 
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triplet pairings is with a relative phase ± J , which leads to 
a novel time-reversal symmetry-breaking pairing stated. 
The investigation of the unconventional Cooper pairing 
symmetry in magnetic dipolar systems was studied by 
the authors 3 -^. We have found that it provides a robust 
mechanism for a novel p- wave (L — 1) spin triplet (S = 1) 
Cooper pairing to the first order in interaction strength. 
It comes directly from the attractive part of the magnetic 
dipolar interaction. In comparison, the triplet Cooper 
pairings in 3 He and solid-state systems come from spin 
fluctuations, which is a second-order effect in interaction 
strengt h 38 i 39 . Furthermore, that pairing symmetry was 
not studied in 3 He systems before in which orbital and 
spin angular momenta of the Cooper pair are entangled 
into the total angular momentum J = 1. In contrast, in 
the 3 He-_B phase 40 , L and S are combined as J = 0, and 
in the 3 He-A phase, L and S are decoupled and J is not 
well-define d 41 i 42 . 

Fermi liquid theory is one of the most important 
paradigms in condensed matter physics on interacting 
fermion o 38 ' 43 . Despite the pioneering papers^^r— , a sys- 
tematic study of the Fermi liquid properties of magnetic 
dipolar fermions is still lacking in the literature. In par- 
ticular, Landau interaction matrices have not been calcu- 
lated, and a systematic analysis of the renormalizations 
from magnetic dipolar interactions to thermodynamic 
quantities has not been performed. Moreover, collective 
excitations in magnetic dipolar ultracold fermions have 
not been studied before. All these are essential parts of 
Fermi liquid theory. The experimental systems of 161 Dy 
and 163 Dy are with a very large hyperfine spin of F = 4p, 
thus the Fermi liquid theory taking into account of all 
the complicated spin structure should be very challeng- 
ing. We take the first step by considering the simplest 
case of spin-i magnetic dipolar fermions which preserve 
the essential features of spin-orbit physics and address 
the above questions. 

In this paper, we systematically investigate the Fermi 
liquid theory of the magnetic dipolar systems includ- 
ing both the thermodynamic properties and the collec- 
tive excitations, focusing on the spin-orbit coupled ef- 
fect. The Landau interaction functions are calculated 
and are diagonalized in the spin-orbit coupled basis. 
Renormalizations for thermodynamic quantities and the 
Pomeranchuk-type Fermi surface instabilities are stud- 
ied. Furthermore, the collective modes are also spin-orbit 
coupled with a topologically non-trivial configuration of 
the spin distribution in momentum space. Their disper- 
sion relation and configurations are analyzed. 

Upon the completion of this paper, we became aware of 
the nice work by Sogo et al£^. Reference constructed 
the Landau interaction matrix for dipolar fermions with 
a general value of spin. The Pomeranchuk instabilities 
were analyzed for the special case of spin i, and collective 
excitations were discussed. Our paper has some overlaps 
on the above topics with Ref. [4J| but with a signifi- 
cant difference, including the physical interpretation of 
the Pomeranchuk instability in the J = 1~ channel and 



our discovery of an exotic propagating spin-orbit sound 
mode. 

The remaining part of this paper is organized as fol- 
lows. The magnetic dipolar interaction is introduced in 
Sec. ITT] The Landau interaction matrix is constructed at 
the Hartree-Fock level and is diagonalized in Sec. IIIII In 
Sec. IIV1 we present the study of the Fermi liquid renor- 
malization to thermodynamic properties from the mag- 
netic dipolar interaction. The leading Pomeranchuk in- 
stabilities are analyzed. In Sec. [V] the spin-orbit coupled 
Boltzmann equation is constructed. We further perform 
the calculation of propagating spin-orbit coupled collec- 
tive modes. We summarize the paper in Sec. IVI1 



II. MAGNETIC DIPOLAR HAMILTONIAN 

We introduce the magnetic dipolar interaction and the 
subtlety of its Fourier transform in this section. 

The magnetic dipolar interaction between two spin-i 
particles located at ri 2 reads 



3(S aa > ■ r)(Sppi ■ f) ,(1) 



where S = |cr; a, a', (3, (3' take values of f and I; f 
r*i — r*2 and f = r/r is the unit vector along r. 
The Fourier transform of Eq. ([TJ is 
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3(S aa > ■ q)(S/3P' ■ q) - S aa > ■ S/3f3< ,(2) 



which depends on the direction along the momentum 
transfer but not its magnitude. It is singular as q — > 0. 
More rigorously, V a p^i a ' (q) should be further multiplied 
by a numeric factor— as 



9(q) = 3 



/ Ji(ge) _ ji(qL) \ 
V qe qL J 



(3) 



where e is a short range scale cut off, and L is the long 
distance cut off at the scale of sample size. The spher- 
ical Bessel function ji(x) shows the asymptotic behav- 
ior ji(x) — > | at x — > 0, and ji{x) i sin(a; — ^) as 
x — > oo. In the long wavelength limit satisfying qe — > 
and qL — > oo, g(q) —> 1 and we recover Eq. @. If q 
is exactly zero, V a p-pi a > = 0, because the dipolar inter- 
action is neither purely repulsive nor attractive, and its 
spatial average is zero. 

The second quantization form for the magnetic dipolar 
interaction is expressed as 



H int = — £ ^ a {k + q)^l(k')V a0 . 0la ,{q) 
k,k' ,q 



(4) 



where V is the volume of the system. The density of 
states of two-component Fermi gases at the Fermi energy 
is A^o = ^rjk , and we define a dimensionless parameter 
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A = No/i 2 . A describes the interaction strength, which 
equals the ratio between the average interaction energy 
and the Fermi energy up to a factor on the order of 1. 



III. SPIN-ORBIT COUPLED LANDAU 
INTERACTION 

In this section, we present the Landau interaction func- 
tions of the magnetic dipolar Fermi liquid, and perform 
the spin-orbit coupled partial wave decomposition. 



A. The Landau interaction function 

Interaction effects in the Fermi liquid theory are cap- 
tured by the Landau interaction function. It describes 
the particle-hole channel forward-scattering amplitudes 
among quasiparticles on the Fermi surface. At the 
Hartree-Fock level, the Landau function is expressed as 



faa',BB>(k,k') = faa'.B0'(q) + /aa' SS'(M') 



(5) 



where k and k' are at the Fermi surface with the mag- 
nitude of kf and q is the small momentum transfer in 
the forward scattering process in the particle-hole chan- 
nel - faa',08'(<l) = v a3,8> a >(q) is the direct Hartree in- 
teraction, and fa a ' B8'(k',k') — —Vaa,a'6'(k — k') is the 
exchange Fock interaction. As q 0, f H is singular, 
thus we need to keep its dependence on the direction of 
q. More explicitly, 



faa',BB'(k\ » ) 



7T/i 



(6) 
(7) 



where the tensor is defined as M aa i t ppi{q) — 3(a aa > ■ 
q){&88' • <?) — 0a a < • uppi and m is the unit vector along 
the direction of the momentum transfer m = . We 



have used the following identity: 



\k-k'\ ' 



3((7 a/ 3' • rh){SBa' ■ m) ~ &a8' ■ VBa' 

3(a aa > ■ rri)(pB8' • rh) - a aa , ■ a BB > 



(8) 



to obtain Eq. ©. 



B. The spin-orbit coupled basis 

Due to the spin-orbit nature of the magnetic dipolar 
interaction, we introduce the spin-orbit coupled partial- 
wave basis for the quasiparticle distribution over the 
Fermi surface following the steps below. 

The 8n aa i (k) is defined as 



where n aa /(k) = {^(k)^ a >{k)) is the Hermitian single- 
particle density matrix with momentum k and satisfies 
n aa i = n* a , a and no(k) is the zero-temperature equilib- 
rium Fermi distribution function ng(k) — 1 — 9{k — kf). 

Sn aa i (fc) is expanded in terms of the particle-hole angular 
momentum basis as 



6n aa ,(k) = ^2Sn Ssz (k)xssz, aa > 



where XSs z ,aa' are the bases for the particle- hole singlet 
(density) channel with 5 = and triplet (spin) channel 
with S = 1, respectively. They are defined as 



X00,aa' — ^aa'i 
XlO, aa' 



(ii) 



which satisfy the orthonormal condition ti(xs s xs's'f) — 

Since quasiparticles are only well defined around the 
Fermi surface, we integrate out the radial direction and 
arrive at the angular distribution, 



5n aa >(k) 



k 2 dk 
J2nf 



Sn aa '(k). 



(12) 



Please note that angular integration is not performed in 
Eq. (fT2|) . We expand 5n aa '{k) in the spin-orbit decou- 
pled bases as 

Sn aa >(k) = ^ 5n LmSS zY Lm (k)xssz,aa>, 

LmSs z 

= 5n *LmS Sz Y Lm(k)Xss a , a a^ ( i3 ) 

LmSs z 

where YL m (k) is the spherical harmonics satisfying the 
normalization condition f dkY£ m (k)YL m (k) = 1. 
We can also define the spin-orbit coupled basis as 

y.}Jz-Ls(k,aa') = y^^(LmSs z \J J z )Y Lm (k)xss z , aa >, 

ms z 

y]j Z ;Ls^aa>) = J2{LmSs z \JJ z )YZ m (k)xl Sz , aa „ 

(14) 

where (LmSs z \JJ z ) is the Clebsch- Gordon coefficient 
and y,jj z -LS satisfies the orthonormal condition of 



6j.J>&LL'8sS>- 



6n aa >(k) = n aa >(k) - S aa >n (k), 



(9) 



(15) 



4 



Using the spin-orbit coupled basis, 8n aa ' (k) is expanded 
as 

5n aa ,(k) = ^2 Srijj^LS yjj z -,Ls(k,ad) 

JJ Z ;LS 

= E 8nl Jt . LS y* JJa . LS (k,acf), (16) 

JJ Z ;LS 

where Sn JJz . LS = J2 mSz (LmSs z \J J z )Sn LmSsz - 

C. Partial- wave decomposition of the Landau 
function 

We are ready to perform the partial-wave decomposi- 
tion for Landau interaction functions. The tensor struc- 
tures in Eqs. ((HJ and only depend on a aa ' and Bppi, 



thus the magnetic dipolar interaction only contributes to 
the spin-channel Landau parameters, i.e., 3=1. In the 
spin-orbit decoupled basis, the Landau functions of the 
Hartree and Fock channels are expanded, respectively, as 

^/fa^3/3'(M') = E Y Lm {k)xisAua') 

Lms z ;L' m' s' z 

(17) 

For later convenience, we have multiplied the density of 
states No and the factor of 1/4tt such that T H - F are di- 
mensionless matrices. Without loss of generality, in the 
Hartree channel, we choose q — z. 

The matrix elements in Eq. (jlT[) are presented below. 
In the Hartree channel, 



T Lmls z ;L'm'W z = "g" ( 2< ^ z ,0 ~ $s z ,±l)$L,0$L',o5m,o8m',o8s z s' z 5 (18) 

and in the Fock channel, 

ttA ( Sll' 5 l+2 ,l' 5l-2,l' 



2 \L(L+1) 3(L+l)(L + 2) 3(L-1)L> 
x / dfl r [6 SzSi - AnY Uz (Q r )Y{ s , (n r )]Y Lm (Q r )Yl, m ,(Sl r ). (19) 



Lmls z ;L'm'le 



The magnetic dipolar interaction is isotropic, thus the spin-orbit coupled basis are the most convenient. In these 
basis, the Landau matrix is diagonal with respect to the total angular momentum J and its z-component J z as 

^/ aa -^(U')= E yjJzMkoca')Fjj zL1 ,jj zLll y\ Jz . vl {k,pp'). (20) 

JJ Z LU 

The matrix kernel Fjj^lx-jj^l'x reads as 

FjJ z Ll-JJ z L'l = -^-5js5 L , 6l>,o(25j^o ~ $J z ,±l) + E ( Lmls ^\ JJ z)( L ' m ' ls 'z\ JJ z) T L m ls z ;L> m >Xs z - ( 21 ) 

ms z :m' s' 



We found that up to a positive numeric factor, the second 
term in Eq. (|21[) is the same as the partial-wave matri- 
ces in the particle-particle pairing channel, which was 
derived for the analysis of the Cooper pairing instability 
in magnetic dipolar systems^. 

However, the above matrix kernel Fjj z l\-,jj z L'\ is 
not diagonal for channels with the same values of J J z 
but different orbital angular momentum indices L and 
L' . Moreover, the conservation of parity requires that 
even and odd values of L do not mix. Consequently, 
Fjj z li-jj z l'i is either diagonalized or reduced into a 
small size of just 2x2. For later convenience of study- 
ing collective modes and thermodynamic instabilities, we 
present below the prominent Landau parameters in some 



low partial- wave channels. Below, we use (J ± J Z LS) to 
represent these channels in which ± represents even and 
odd parities, respectively. 

The parity odd channel of J = 0~ only has one possi- 
bility of ((Toil) in which 

7T 

^0-011;0-011 — 

There is another even parity density channel with J = 
+ , i.e., (0 + 000), which receives contribution from short 
range s-wave interaction but no contribution from the 
magnetic dipolar interaction at the Hartree-Fock level. 
The parity odd channel of J = 1~ only comes from 
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(1 J z ll) in which 



(23) 



Another channel of J = 1~, i.e., (1~ J z 10), channel from 
the p-wave channel density interactions, which again re- 
ceives no contribution from magnetic dipolar interaction 
at the Hartree-Fock level. These two J = 1~ modes are 
spin- and charge-current modes, respectively, and thus, 
do not mix due to their opposite symmetry properties 
under time-reversal transformation. 

We next consider the even parity channels. The 
J = 1 + channels include two possibilities of (JJ Z LS) = 
(1 + J 2 01), (1 + J z 21). The former is the ferromagnetism 
channel, and the latter is denoted as the ferronematic 
channel in Refs. @ and (24[. Due to the spin-orbit na- 
ture of the magnetic dipolar interaction, these two chan- 
nels are no longer independent but are coupled to each 
other. Because the Hartree term breaks the rotational 
symmetry, the hybridization matrices for J z = 0, ±1 are 
different. For the case of J z — 0, it is 



1+0 



-Fi001;1001 -Fl001;1021 
-^1021:1001 ^1021:1021 



ttA / 8 V2\ 
12 \V2 1 J ' 

(24) 



whose two eigenvalues and their associated eigenvectors 
are 

w\ +0 = 0.69ttA, Vi +0 = (0.98,0.19) T , 

w \ +0 = 0.06vrA, ?/> 2 1+0 = (-0.19,0.98) T . (25) 

The hybridization is small. For the case of J z = ±1, the 
Landau matrices are the same as 



F 1+1 = 



_ I i*1101;1101 -Fll01;1121 \ _ ^ ( ~ 4 \[2 



21:1101 



1 Fll21;1121 J 12 \ V2 1 



(26) 



Again the hybridization is small as shown in the eigen- 
values and their associated eigenvectors 



1+1 



-0.37ttA, = (0.97,-0. 25)^ 



w. 



1+1 



0.12ttA, 



(0.25, 0.97) 5 



(27) 



Landau parameters, or, matrices, in other high partial- 
wave channels are neglected, because their magnitudes 
are significantly smaller than those above. 

We need to be cautious on using Eqs. and (|21)|) in 
which the Hartree contribution of Eq. [5] is taken. How- 
ever, Eq. (|6|) is valid in the limit q <C kt but should 
be much larger than the inverse of sample size 1/L. It 
is valid to use Eqs. (|2"4l and (|2"6l when studying the 
collective spin excitations in Sec. [V] below. However, 
when studying thermodynamic properties, say, magnetic 
susceptibility, under the external magnetic-field uniform 
at the scale of L, the induced magnetization is also uni- 
form. In this case, the Hartree contribution is suppressed 



to zero, thus the Landau matrices in the J — 1 + channel 
are the same for all the values of J z as 



FlJ x 01;lJ z 01 -FlJ.01;U^21 
FlJ z 21;lJ z 01 -Fl J. 21; 1.4 21 



than 



ttA 
12 



V2 
V2 1 



(28) 



In this case, the hybridization between these two chan- 
nels is quite significant. The two eigenvalues and their 
associated eigenvectors are 



(29) 



IV. THERMODYNAMIC QUANTITIES 

In this section, we study the renormalizations for ther- 
modynamic properties by the magnetic dipolar interac- 
tion and investigate the Pomeranchuk-type Fermi surface 
instabilities. 



A. Thermodynamics susceptibilities 

The change in the ground-state energy with respect 
to the variation in the Fermi distribution density matrix 
include the kinetic and interaction parts as 



SE SE km , SEi, 



V V V 



(30) 



The kinetic-energy variation is expressed in terms of the 
angular distribution of 5n aa i (k) as 



8E kin Air 



V 



J dkSn aa '{k)Sn a ' a (k) 



TF Y] §n lmSsJ n LmSs z , (31) 
iV ' — 



LmSSz 



where the units of Sns Sz (k) and 8nLmSs z are the same 
as the inverse of the volume. The variation in the inter- 
action energy is 

SE y nt = ^ X! J J dkdk' y aa >,00'{k,k')6n a i a (k)6np'p(k') 

= 2 ^ Sn L m SsjLTnSs z ,L'm'Ss'Jn* L , m , s ^. 

Lms z L'm' s' z ;S 

(32) 

Adding them together and changing to the spin-orbit 
coupled basis, we arrive at 

~V = W 51 Sn *JJ z ;LS M JJ z LS;JJ z L>S$njj z ; L ,s, (33) 
JJ Z ;LL'-S 
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where the matrix elements are 

Mjj z LS;JJ z L'S = 5lL> + FjJ z LS;JJ z L'S- (34) 

In the presence of the external field Hjj^ls, the ground 
state energy becomes 

~V = 16n {j~ ^ Snjj x . LS Mjj zLS ;jj t L>sSnjj g . L > s 

X ° JJ Z LL'S 

- ^2 h J-hLsSnj Jz . Ls y (35) 

JJ Z LS 

where Xo = N is the Fermi liquid density of states. At 
the Hartree-Fock level, Nq receives no renormalization 
from the magnetic dipolar interaction. The expectation 
value of Snjj z i J s is calculated as 

SnjJ^LS = XQ J2( M }jJzLS;JJzL'S h JJ.L'S- (36) 
V 

For the J = 1 + channel, M^ 1 m I — Fi+ thmW U P to 
first order of A in the case of A <C 1. As a result, the exter- 
nal magnetic field h along the z axis not only induces the 
^-component spin polarization, but also induces a spin- 
nematic order in the channel of (J + J Z LS) = (1 + 021), 
which is an effective spin-orbit coupling term as 

SH = ^7r\h^Uk){[(k 2 -3k 2 z )a z 

k 

- 3k z (k x a- x + ky(jy)] \ipp(k). (37) 

Apparently, this term breaks time-reversal symmetry, 
and thus cannot be induced by the relativistic spin-orbit 
coupling in solid states. This magnetic field induced spin- 
orbit coupling in magnetic dipolar systems was studied 
by Fregoso et al^^ 

B. Pomeranchuk instabilities 

Even in the absence of external fields, Fermi surfaces 
can be distorted spontaneously known as Pomeranchuk 
instabilities 45 . Intuitively, we can imagine the Fermi sur- 
face as the elastic membrane in momentum space. The 
instabilities occur if the surface tension in any of its 
partial-wave channels becomes negative. In the magnetic 
dipolar Fermi liquid, the thermodynamic stability condi- 
tion is equivalent to the fact that all the eigenvalues of 
the matrix M JJzLS -jj z l'S are positive. 

We next check the negative eigenvalues of the Landau 
matrix in each partial-wave channel. Due to the absence 
of external fields, the Pomeranchuk instabilities are al- 
lowed to occur as a density wave state with a long wave 
length q — > 0. For the case of J = f + , it is clear that in 
the channel of J z = ±1, the eigenvalue w\ +1 in Eq. ([27)1 
is negative and the largest among all the channels. Thus 
the leading channel instability is in the (JJ Z ) = (1 + ± 1) 



channel, which occurs at w\ 1 < —1, or, equivalently, 
A > A^ +1 = 0.86. The corresponding eigenvector shows 
that it is mostly a ferromagnetism order parameter with 
small hybridization with the ferronematic channel. A re- 
pulsive short-range s wave scattering, which we neglected 
above will enhance ferromagnetism and, thus, will drive 
AJ +1 to a smaller value. The wavevector q of the spin 
polarization should be on the order of 1/L to minimize 
the energy cost of twisting spin, thus, essentially exhibit- 
ing a domain structure. The spatial configuration of the 
spin distribution should be complicated by actual bound- 
ary conditions. In particular, the three- vector nature of 
spins implies the rich configurations of spin textures. An 
interesting result is that the external magnetic field actu- 
ally weakens the ferromagnetism instability. If the spin 
polarization is aligned by the external field, the Landau 
interaction matrix changes to Eq. (l28l) . The magnitude 
of the negative eigenvalue is significantly smaller than 
that of Eq. (j!?6"|) . As a result, an infinitesimal external 
field cannot align the spin polarization to be uniform but 
a finite amplitude is needed. 

For simplicity, we only consider ferromagnetism with 
a single plane wave vector q along the z axis, then the 
spin polarization spirals in the xy-plane. Since q ~ 1/L, 
we can still treat a uniform spin polarization over a dis- 
tance large comparable to the microscopic length scale. 
Without loss of generality, we set the spin polariza- 
tion along the x axis. As shown in Ref. [24|, ferro- 
magnetism induces ferronematic ordering. The induced 
ferronematic ordering is also along the x axis, whose 
spin-orbit coupling can be obtained based on Eq. (|37[) 
by a permutation among components of k as H' so (k) oc 
(fc 2 -3k x )a x -3k x (k y a y + k z a z ). According to Eq. (p7|i . 
ferromagnetism and ferronematic orders are not strongly 
hybridized, the energy scale of the ferronematic SO cou- 
pling is about 1 order smaller than that of ferromag- 
netism. An interesting point of this ferromagnetism is 
that it distorts the spherical shape of the Fermi surface as 
pointed by Fregoso and Fradkir*21. This anisotropy will 
also affect the propagation of Goldstone modes. Further- 
more, spin waves couple to the oscillation of the shape of 
Fermi surfaces bringing Landau damping to spin waves. 
This may result in non-Fermi liquid behavior for fermion 
excitations, and will be studied in a later paper. This ef- 
fect in the nematic symmetry-breaking Fermi liquid state 
has been extensively studied before in the literature 4 ^—. 

The next subleading instability is in the J = 1~ chan- 
nel with L = 1 and S = 1 as shown in Eq. (1231 . which 
is a spin-current channel. The generated order parame- 
ters are spin-orbit coupled. For the channel of J z = 0, 
the generated SO coupling at the single-particle level ex- 
hibits the three-dimensional (3D) Rashba type as 

H so ,i- = \ n z\^2 / ipl l {k){k x ay - kya x ) a0 ip^(k), (38) 

where \n z \ is the magnitude of the order parameter. The 
same result was also obtained recently in Ref. 13 In 
the absence of spin-orbit coupling, the L = S = 1 chan- 
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nel Pomeranchuk instability was studied in Refs. |52j | 
and [531 ] . which exhibits the unconventional magnetism 
with both isotropic and anisotropic versions. They 
are particle-hole channel analogies of the p-wave triplet 
Cooper pairings of 3 He isotropic B and anisotropic A 
phases, respectively. In the isotropic unconventional 
magnetic state, the total angular momentum of the or- 
der parameter is J = 0, which exhibits the k ■ cr-type 
spin-orbit coupling. This spin-orbit coupling is gener- 
ated from interactions through a phase transition and, 
thus, was denoted as the spontaneous generation of spin- 
orbit coupling. In Eq. (|3"5f . the spin-orbit coupling that 
appears at the mean-field single-particle level cannot be 
denoted as spontaneous because the magnetic dipolar in- 
teraction possesses the spin-orbit nature. Interestingly, 
in the particle-particle channel, the dominant Cooper 
pairing channel has the same partial-wave property of 
L = S = J = lZL. 

The instability in the J = 1~ (spin current) channel 
is weaker than that in the 1 + (ferromagnetism) channel 
because the magnitude of Landau parameters is larger 
in the former case. The 1~ channel instability should 
occur after the appearance of ferromagnetism. Since 
spin-current instability breaks parity, whereas, ferromag- 
netism does not, this transition is a genuine phase tran- 
sition. For simplicity, we consider applying an external 
magnetic field along the z axis in the ferromagnetic state 
to remove the spin texture structure. Even though the 
J = 1 + and channels share the same property un- 



der rotation transformation, they do not couple at the 
quadratic level because of their different parity proper- 
ties. The leading-order coupling occurs at the quartic 
order as 



SF = /3 x (n -n)(S- S) + (3 2 \n x 5| 2 



(39) 



where ft and S represent the order parameters in the 
J = 1~ and 1 + channels, respectively. fi\ needs to be 
positive to keep the system stable. The sign of $1 de- 
termines the relative orientation between n and S. It 
cannot be determined purely from the symmetry analy- 
sis but depends on microscopic energetics. If (3 2 > 0, it 
favors n II S, and n _L S is favored at /3 2 < 0. 



V. THE SPIN-ORBIT COUPLED COLLECTIVE 
MODES 

In this section, we investigate another important fea- 
ture of the Fermi liquid, the collective modes, which again 
exhibit the spin-orbit coupled nature. 



A. Spin-orbit coupled Boltzmann equation 

We employ the Boltzmann equation to investigate the 
collective modes in the Fermi liquid stated 



J 



9 n(f, k, t) - i [e(r, k, t) , n(f, k, + 



ot 



dh 



dn 



de(f,k,t) dn(f,k,t) 



dri 



dh J 



0, (40) 



where n aa i (r, fc, t) and e aa ' (r, k, t) are the density and en- we arrive at 
ergy matrices for the coordinate (r, k) in the phase space 
and [, ] and {, } mean the commutator and anticommuta- Sn aa ' (k) - 
tor, respectively. Under small variations in n aa >(f,k,t) 
and e aa i (r, k,t), 



1 cos ( 



2 s — cos ( 



x 8n pf}l (k') = 0, 



k'^faa',f}fl'{k 7 k') 



(43) 



n aa >(r,k,t) = n {k)5 aa ' + Sn aa '(r,k,t), 

f d 3 k' - - 

e„„/(r,fc,t) = e(fc)<W + / 7^35 faa',/30'(k, k') 



(2tt) 2 



x 8npi3>{k') 



(41) 



the above Boltzmann equation can be linearized. Plug- 
ging the plane-wave solution of 



bn aa , (f, fc, t) = 8n aa , (k)e t{q ' r ~ 



LOt) 



(42) 



where s is the dimensionless parameter uj/(vfq). The 
propagation direction of the wavevector q is defined along 
the z-direction. 

In the spin-orbit decoupled basis defined as 8riLrnSs z 
in Sec. 1111 B| the linearized Boltzmann equation becomes 

8riLmSs z + ^LL';m{s)FL'm'Ss z ;L"m"Ss'j8nL"m"Ss'J = 0, 

' (44) 

where Q,ll'(s) is equivalent to the particle-hole channel 
Fermi bubble in the diagrammatic method as 



LL':m 



^^WW*)^^. (45) 
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For later convenience, we present Q.LL';m m several chan- 
nels of LL' and m as follows 

£Vo(s) = 1- Jln|— ti|+i|se(s<l), 
2 1 — s 2 



^10;o(s) = ^01:0 = \/3s^00;0 («) , 

^ii;o(s) = 1 + 3s 2 fi o;o(s), 



^ll;l(s) = ^ll;-l(s) 



1 r 

2 



1 - 3(1 - s 2 )n 00 . {s) 



(46) 



Equation (|44j) can be further simplified by using the 
spin-orbit coupled basis Snjj Zi i,s defined in Sec. 1111 Bl 

3njJz-,LS + ^ Kjj xLS -j>JzL's{s)Fj/j xL /s-j>JzL"S 

J';LL< 

x 6nj>j aL »s=0, (47) 
where the matrix kernel Kjj^ls;J>j z l'S reads 
K.j.j z ls;J>JzL>s{s) = y^^(LmSs z \JJ z )(L'mSs z \J'J z ) 

ms z 

x n LL ,, m (s). (48) 

B. The spin-orbit coupled sound modes 

Propagating collective modes exist if Landau parame- 
ters are positive. In these collective modes, interactions 
among quasiparticles rather than the hydrodynamic col- 
lisions provide the restoring force. Because only the spin 
channel receives renormalization from the magnetic dipo- 
lar interaction, we only consider spin channel collective 



modes. The largest Landau parameter is in the (1 + 001) 
channel in which the spin oscillates along the direction 
of q. The mode in this channel is the longitudinal spin 
zero sound. On the other hand, due to the spin-orbit 
coupled nature, the Landau parameters are negative in 
the transverse spin channels of (1 + ±1 ±1), and thus 
no propagating collective modes exist in these channels. 
The hybridization between (1+001) and (1 + 021) is small 
as shown in Eq. (1251) . and the Landau parameter in the 
(1 + 021) channel is small, thus, this channel also is ne- 
glected below for simplicity. 

Because the propagation wave vector q breaks the par- 
ity and 3D rotation symmetries, the (1 + 001) channel cou- 
ples to other channels with the same J z . As shown in 
Eq. (|47l) . the coupling strengths depend on the magni- 
tudes of Landau parameters. We truncate Eq. (14"T1) by 
keeping the orbital partial-wave channels of L = and 
L = 1 because Landau parameters with orbital-partial 
waves L > 2 are negligible. There are three channels with 
L = S = 1 as (0-011), (1^011), and (2-011). We further 
check the symmetry properties of these four modes un- 
der the reflection with respect to any plane containing q. 
The mode of (1 — 011) is even and the other three are odd, 
thus it does not mix with them. The Landau parameter 
in the (2-011) channel is calculated as ^A, which is 1 
order smaller than those in (1+001) and (1-001), thus 
this channel is also neglected. We only keep these two 
coupled channels (1+001) and (1~001) in the study of 
collective spin excitations. 

The solution of the two coupled modes reduces to a 
2x2 matrix linear equation as 



( 1 + fioo ; o(s)-Fl001;1001 s ^OO;o( s )-Pboil;0011 \ f Sn W Qi \ ^ 
\ s ^OO;o( s )^1001;1001 1 + ^00;o( s )-Po011;0011 / \ ^0011 / 

where the following relations are used 

A"iooi ; iooi(s) = f2oo;o(s) 

Xiooi;oon(«) = #ooii;iooi («) = (0010|10)<1010|00)Q O i ; o(s) = sfioojoOO 

#ooii,oon(s) = ^|(l™l-m|00)| 2 Oii ;m (s) = -Oii ; o(s) + -0 U; i(s) =n (s0 , {s). (50) 

771 
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The condition of the existence of nonzero solutions of Eq. 
(|49)l becomes 

(1 - .s 2 )^ 0;0 ( S ) + 2n 00 ;o(s& + 4- = 0, (51) 

w here F + = {F 10Q 1 . 1001 + F 00 ii;ooii)/2 and F x = 

\f -Fio01:1001^0011;0011- 



Let us discuss several important analytical properties 
of its solutions. In order for collective modes to prop- 
agate in Fermi liquids, its sound velocity must satisfy 
s > 1, otherwise it enters the particle- hole continuum 
and is damped, a mechanism called Landau damping. 



FIG. 1: (Color online) The sound velocity s in the unit of 
Vf v.s. the dipolar coupling strength A. At < A < 1, 
s(A) ~ 1 + + . On the order of A S> 1, s(A) becomes linear 
with the slope indicated in Eq. (|56|) . 



We can solve Eq. (|5Tj) as 



00;0V 



(* 2 



I) Ft 



(52) 



Only the expression of the S1q . (s) is consistent with s > 
1 and is kept. The other branch has no solution of the 
propagating collective modes. 

Let us analytically check two limits with large and 
small values of A, respectively. In the case of < A -C 1 
such that s — > 1 + 0+ , Eq. flSTJ reduces to 



^oo ; o(sa<i) ~ 1 



1 



In 2 



1 



Its sound velocity solution is 



S A<£1 



l + 2e 



ln(s - 1) 



= 1 + 2e~ 



1 



2F+ 



(53) 



(54) 



The eigenvector can be easily obtained as ^^(1, 
which is an equal mixing between these two modes. On 
the other hand, in the case of A 3> 1, we also expect 
s>l, and thus Eq. (fSTj) reduces to 



^00;o(sA»l) ~ — — TT-n 



sF x 



1 

3s 2 ' 



whose solution becomes 



S A»1 



3 



3^3 



(55) 



(56) 



In our case, Fiooi is larger than Fqqh but is on 
the same order. The eigenvector can be solved as 
i \J TtooT) 



-(VF, 



0011, 



in which the weight of the 



V2F+ ' 

(0011) channel is larger. 

The dispersion of the sound velocity s with respect to 
the dipolar interaction strength A is solved numerically 
as presented in Fig. [TJ Collective sound excitations exist 
for all the interaction strengths with s > 1. In both 
limits of <C A <C 1 and A 3> 1, the numerical solutions 



k./k„ 



FIG. 2: (Color online) The spin configuration [Eq. ((57)) ] of 
the zero-sound mode over the Fermi surface shows hedgehog- 
type topology at A = 10. The common sign of u\ and U2 is 
chosen to be positive, which gives rise to the Pontryagin index 
+ 1. Although the hedgehog configuration is distorted in the 
z component, its topology does not change for any values of 
A describing the interaction strength. 



agree with the above asymptotic analysis of Eqs. (|54j) 
and (|56"|) . In fact, the linear behavior of s(A) already 
appears at A ~ 1, and the slope is around 0.6. For all the 
interaction strengths, the (1 + 001) and (0~011) modes 
are strongly hybridized. 

This mode is an oscillation of spin-orbit coupled Fermi 
surface distortions. The configuration of the (0~011) 
mode exhibits an oscillating spin-orbit coupling of the 
k ■ a type. This is the counterpart of the isotropic un- 
conventional magnetism, which spontaneously generates 
the k ■ a- type couplin g 52 ! 53 . The difference is that, here, 
it is a collective excitation rather than an instability. It 
strongly hybridizes with the longitudinal spin mode. The 
spin configuration over the Fermi surface can be repre- 
sented as 



s(r, k, t) 



u 2 sin Oj: cos 0£ 
U2 sin 9^ sin (f>^ 
u 2 cos 0g + u\ 



J(q-f-sqvft) 



(57) 



where (ui, U2) T is the eigenvector for the collective mode. 
We have checked that for all the values of A, \u 2 \ > \ui\ 
is satisfied with no change in their relative sign, thus the 
spin configuration as shown in Fig. [5] is topologically 
non-trivial with the Pontryagin index ±1 which periodi- 
cally flips the sign with time and the spatial coordinate 
along the propagating direction. It can be considered as 
a topological zero sound. 



VI. CONCLUSIONS 

To summairze, we have presented a systematic study 
on the Fermi liquid theory with the magnetic dipolar 
interaction, emphasizing its intrinsic spin-orbit coupled 
nature. Although this spin-orbit coupling does not ex- 
hibit at the single-particle level, it manifests in various 
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interaction properties. The Landau interaction function 
is calculated at the Hartree-Fock level and is diagonal- 
ized by the total angular momentum and parity quan- 
tum numbers. The Pomeranchuk instabilities occur at 
the strong magnetic dipolar interaction strength gener- 
ating effective spin-orbit coupling in the single-particle 
spectrum. 

We have also investigated novel collective excitations 
in the magnetic dipolar Fermi liquid theory. The Boltz- 
mann transport equations are decoupled in the spin-orbit 
coupled channels. We have found an exotic collective ex- 



citation, which exhibits spin-orbit coupled Fermi surface 
oscillations with a topologically nontrivial spin configu- 
ration, which can be considered as a topological zero- 
sound-like mode. 



Acknowledgments 

Y. L. and C. W. were supported by the AFOSR YIP 
program and NSF-DMR-1105945. 



1 S. Ospelkaus, K. K. Ni, M. H. G. de Miranda, B. Neyen- 
huis, D. Wang, S. Kotochigova, P. S. Julienne, D. S. Jin, 
and J. Ye, Faraday Discuss. 142, 351 (2009).. 

2 K. K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, 
B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, 
D. S. Jin, and J. Ye, Science 322, 231 (2008). 

3 A. Chotia, B. Neyenhuis, S. A. Moses, B. Yan, J. P. Covey, 
M. Foss-Feig, A. M. Rey, D. S. Jin, and J. Ye, Phys. Rev. 
Lett. 108, 080405 (2012). 

4 T. Sogo, L. He, T. Miyakawa, S. Yi, H. Lu, and H. Pu, 
New J. Phys. 11, 055017 (2009). 

5 T. Miyakawa, T. Sogo, and H. Pu, Phys. Rev. A 77, 061603 
(2008). 

6 B. M. Fregoso, K. Sun, E. Fradkin, and B. L. Lev, New J. 
Phys. 11, 103003 (2009). 

7 C. K. Chan, C. Wu, W. C. Lee, S. Das Sarma, Phys. Rev. 
A 81, 023602 (2010), 

8 S. Ronen and J. L. Bohn, Phys. Rev. A 81, 033601 (2010). 

9 C. Lin, E. Zhao, and W. V. Liu, Phys. Rev. B 81, 045115 
(2010). 

10 Y. Yamaguchi, T. Sogo, T. Ito, T. Miyakawa, Phys. Rev. 
A 82, 013643 (2010). 

11 Q. Li, E. H. Hwang, S. Das Sarma, Phys. Rev. B 82, 
235126 (2010). 

12 C. M. Chang, W. C. Shen, C. Y. Lai, P. Chen, D. W. 
Wang, Phys. Rev. A 79, 053630 (2009). 

13 J. N. Zhang, S. Yi, Phys. Rev. A 81, 033617 (2010); J. N. 
Zhang, S. Yi, Phys. Rev. A 80, 053614 (2009). 

14 K. Sun, C. Wu, S. Das Sarma, Phys. Rev. B 82, 075105 
(2010). 

15 S. H. Youn, M. Lu, U. Ray, and B. L. Lev, Phys. Rev. A 
82, 043425 (2010). 

16 M. Lu, S. H. Youn, and B. L. Lev, Phys. Rev. Lett. 104, 
063001 (2010). 

17 M. Lu, N. Q. Burdick, S. H. Youn, B. L. Lev, Phys. Rev. 
Lett. 107, 190401 (2011). 

18 T. Koch, T. Lahaye, J. Metz, B. Frohlich, A. Griesmaier, 
and T. Pfau, Nature Physics 4, 218 (2008). 

19 T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. 
Pfau, Rep. Prog. Phys. 72, 126401 (2009). 

20 T. Lahaye, J. Metz, T. Koch, B. Frohlich, A. Griesmaier, 
and T. Pfau, Pushing the Frontiers of Atomic Physics, 1, 
160 (2009). 

21 C. Menotti, M. Lewenstein, T. Lahaye, and T. Pfau, in 
Dipolar Interaction in Ultra-Cold Atomic Gases, edited by 
A. Campa, A. Giansanti, G. Morigi and F. S. Labini, AIP 
Conf. Proc. No. 970 (AIP, New York, 2008), p. 332. 



22 J. Quintanilla, S. T. Carr, and J. J. Betouras, Phys. Rev. 
A 79, 031601 (R) (2009). 

23 B. M. Fregoso, E. Fradkin, Phys. Rev. Lett. 103, 205301 
(2009) 

24 B. M. Fregoso, E. Fradkin, Phys. Rev. B 81, 214443 (2010). 

25 L. You and M. Marinescu, Phys. Rev. A 60, 2324 (1999). 

26 M. A. Baranov, M. S. Mar'enko, V. S. Rychkov, and G. V. 
Shlyapnikov, Phys. Rev. A 66, 013606 (2002). 

27 M. A. Baranov, L. Dobrek, and M. Lewenstein, Phys. Rev. 
Lett. 92, 250403 (2004). 

28 M. A. Baranov, Physics Reports 464, 71 (2008). 

29 G. M. Bruun, E. Taylor, Phys. Rev. Lett. 101, 245301 
(2008). 

30 J. Levinsen, N. R. Cooper, and G. V. Shlyapnikov, Phys. 
Rev. A 84, 013603 (2011). 

31 A. C. Potter, E. Berg, D. W. Wang, B. I. Halperin, and E. 
Demler, Phys. Rev. Lett. 105, 220406 (2010). 

32 R. M. Lutchyn, E. Rossi, and S. Das Sarma, Phys. Rev. A 
82, 061604 (2010). 

33 C. Wu, and J. E. Hirsch, Phys. Rev. B 81, 020508 (R) 

(2010) . 

34 K. V. Samokhin, and M. S. Mar'enko, Phys. Rev. Lett. 97, 
197003 (2006). 

35 T. Shi, J. N. Zhang. C. P. Sun and S. Yi, Phys. Rev. A 82, 
033623, (2010). 

36 B. Kain, Hong Y. Ling, Phys. Rev. A 83, 061603(R), 

(2011) . 

37 Y. Li, C. Wu, Scientific Reports 2, 392 (2012). 

38 A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975). 

39 G. E. Volovik, "The Universe in a Helium Droplet" , (Ox- 
ford University Press, Oxford 2009). 

40 R. Balian, N. R. Werthamer, Phys. Rev. 131, 1553 (1963). 

41 P. W. Anderson, P. Morel, Phys. Rev. 123, 1911 (1961). 

42 W. F. Brinkman, J. W. Serene, and P. W. Anderson, Phys. 
Rev. A 10, 2386 (1974). 

43 W. Negele and H. Orland, Quantum Many-Particle Sys- 
tems, (Perseus Books, New York, 1988). 

44 T. Sogo, M. Urban, P. Schuck, T. Miyakawa, Phys. Rev. 
A 85, 031601(R) (2012). 

45 I. J. Pomeranchuk, Sov. Phys. JETP 8, 361 (1959). 

46 V. Oganesyan, S. A. Kivelson, and E. Fradkin, Phys. Rev. 
B 64, 195109 (2001). 

47 M. Garst and A. V. Chubukov, Phys. Rev. B 81, 235105 
(2010). 

48 M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 075127 
(2010) 

49 P. Woffle, and A. Rosch, Journal of Low Temperature 



11 



Physics, 147, 165 (2007). 

C. A. Lamas, D. C. Cabra, and N. Grandi, Phys. Rev. B 
80, 075108 (2009). 

C. A. Lamas, D. C. Cabra, and N. E. Grandi, Int. J. of 
Mod. Phys. B 25, 3539 (2011). 



° 2 C. Wu, and S. C. Zhang, Phys. Rev. Lett. 93, 
036403(2004). 

53 C. Wu, K. Sun, E. Fradkin, and S. C. Zhang, Phys. Rev. 
B 75, 115103 (2007). 



